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GSH 006—15+7: A local Galactic supershell featuring 
transition from Hi emission to absorption 
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ABSTRACT 

We report on the discovery of a new Galactic supershell, GSH 006—15+7, from the Galactic 
All Sky Survey data. Observed and derived properties are presented and we find that 
GSH 006—15+7 is one of the nearest physically large supershells known, with dimensions 
of ~ 780 x 520 pc at a distance of ~ 1.5 kpc. The shell wall appears in Hi emission at 
b < —6.5° and in Hi self-absorption (HlSA) at b > —6.5°. We use this feature along with 
HlSA diagnostics to estimate an optical depth of r ~ 3, a spin temperature of ~ 40 K and 
a swept-up mass of M ~ 3 x 10 6 M . We also investigate the origin of GSH 006—15+7, 
assessing the energy contribution of candidate powering sources and finding evidence in 
favour of a formation energy of ~ 10 52 ergs. We find that this structure provides evidence 
for the transfer of mass and energy from the Galactic disk into the halo. 

Key words: ISM: bubbles - ISM: structure - Galaxy: structure - radio lines: ISM - 
radiative transfer - stars: winds 



1 INTRODUCTION 

The pervasiveness of neutral hydrogen (Hi) in the universe 
has ensured its centrality in both Galactic and extragalactic 
studies, from star-forming regions to supernovae to disc-halo 
interaction. La rge shell-like structures detected in Hi, known 
as supershells (|Heileslll976r h are an important driving force in 
the Galactic ecosystem, circulating material within the disc, 
evacuating dense regions and encouraging n e w star formation 
(ICox fe Sm ith 1974; lMcKee fe Os triker 1973 ;lBruhweiler et al.1 
Il98d : iMcCrav fe Kafato s 1987; Weis z et al.ll2009h . It is theo- 
rised that supershells also provide a pivotal link between th e 
disc of a galaxy and the material in its halo (|Dove et al.ll2000[ ), 
indicating that there should be shared properties between halo 
clouds that are Galactic in origin a nd supershells on the verge 
of blo wing out of the Galactic disc (|Lockmanll2002l ; I Ford et all 
l2008h . 

Supershells are believed to be formed through the com- 
bined effects of multiple supernovae exploding in a region 
on the order of parsecs to hun dreds of pa rsecs in size, of- 
ten around an OB association (|Heileslll979h . It is not clear 
whether this combined effort occurs during one generation of 
high-mass stars or whether it is t he outcome of m any gen- 
erations of stars within a region (|Oev et al.1 l2005h , but the 
resulting supershell of expanding neutral hydrogen undoubt- 
edly requires significant energy input, from 10 52 to as much 



as 10 54 ergs, in order to become a physical structure hundreds 
of parsecs in size. The lifetime of a typical supershell is of the 
order of 10 7 years (| Tenor io-Tagle fe Bodenheimerlll988h , mean- 
ing that these structures signifi cantly outlast their paren t su- 
pernova remnants (~ 10 5 vear sJGuseinov fe Yu sifov 198 and 
can t r ace past star format ion (| Weisz et al. 1 20091 : lAnathpindikal 
l201ll : ICannon et al.l l201lh as well as induce new star forma- 
tion (|Ortega et alJ l2Q09). It is generally difficult to identify the 
stars which formed an old supershell, given that most of the 
high-mass stars will have been e xtinguished as supernova e in 
the process of forming the shell (|McCrav fe Kafatod [l987h . If 
a potential powering source is identified, it is possible to assess 
whether the formation of the vacuou s cavity can be accounted 
for by its winds and/or supernova e ( Bro wn et al.1 Il995l : IShullI 
ll995l : lMcClure-Griffiths et al. II 200 ltlSmithl 20061 ): in most cases, 
however, it is impossible to find a source region in an old shell 
(t > 5 Myr). 

We report on the discovery of a Galactic supershell, 
GSH 006—15+7, an object which exhibits a coherent struc- 
ture over an angular size of ~ 25°. In Section [21 we describe 
its observed and derived properties in the context and theory 
of known Galactic supershells and investigate mult i wavelength 
data in order to characterise GSH 006—15+7. In Section [3l 
we use the unusual property of transition from Hi emission to 
self-absorption visible in the walls of the shell to derive spin- 
temperature, kinetic temperature and optical depth, which al- 
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GALACTIC Longitude 

Figure 1. Hi image of GSH 006—15+7 at the systemic velocity of 7 
km s _1 . The brightness ranges from 10—80 K on a linear scale. Each 
velocity slice has a width of ~ 1 km s _1 . 



lows us to better estimate the mass of the shell. Finally, in 
Section [4] we look for potential powering sources that could 
account for the energy required to form the shell. This super- 
shell presents an interesting case study as one of the nearest 
large supershells. It highlights the potential dangers of assum- 
ing negligible optical depth when estimating mass from column 
density as well as being a possible example of a supershell in 
the transition stage between expansion and blowout. Based on 
the candidate powering sources, we also find evidence in favour 
of an origin for the supershell in the Sagittarius spiral arm. 



2 GSH 006-15+7: MORPHOLOGY AND 
PHYSICAL PROPERTIES 

GSH 006—15+7 was discovered within the Galactic All-Sky 
Sur vey (GASS) data set, a recent southern-sky Hi line d ata sur- 
vey (|McClure-Griffiths et al.ll2009l : lKalberla et al.ll2010h . GASS 
covers the entire sky with declination 5 ^ 1° at an angular res- 
olution of no 16 arcmin, a velocity resolution of ^ 1 km s _1 
over the velocity range —400 ^ v ^ 500 km s _1 and with an 
RMS brightness temperature of ~ 55 mK. The data set was 
produced using the 21 cm multibeam receiver on the Parkes 
64 m radio telescope. 



2.1 General morphology 

The shell is named GSH 006—15+7 after its central longi- 
tude (no 6°), central latitude (~ —15°) and systemic veloc- 
it£l (vlsr - 7 km s" 1 ). GSH 006-15+7 is an example of 
a well-defined supershell of large angular size, comparable to 
those described in previous surveys for Galactic supershells 
(IHeiled 1 19791 : IHuI Il98ll : IHeiled 1 19841 : iMcClure-Griffiths et al.1 
I2002T ). It extends over longitudes of / = [356,16]° and lati- 
tudes of b — [—28,2]°, and is evident over the velocity range 
i^lsr ~ 5 — 11 km s _1 . Fig. [1] shows the shell in the ve- 
locity slice corresponding to its assumed systemic velocit}0, 
i^lsr ~ 7 km s _1 , while Fig. [2] shows its morphology over the 
multiple velocity slices in which it is apparent. We expect that 
there is emission due to the shell at velocities below 5 km s _1 
which cannot be discerned due to the confusion of local emis- 
sion. Table [1] provides a summary of all observed and derived 
properties. 

The shell is teardrop-shaped with its apex crossing 
through the Galactic plane and extends to a latitude of —25°. 
This teardrop morphology suggests an origin in the Galac- 
tic plane, with material elongated in the latitude direction. 
GSH 006—15+7 is an asymmetric structure as there is no con- 
vincing evidence of an upper-latitude counterpart with respect 
to the Galactic plane. The region the shell occupies is con- 
fused due to its low velocity, resulting in considerable overlap 
with local emission, but the shell remains quite well-defined 
over a large angular size. We examine the structure of the shell 
at high latitude for evidence of fragmentation and breakout 
(Section 12. 4p , and find an apparent break in the shell wall at 
/ no 10°, b no —25°, particularly at the lowest velocities (see Fig. 
[3]). There is also a transition from Hi emission to absorption 
on the upper right wall of the shell towards the Galactic plane 
over b = [—7.0,-3.0]°, typically indicative of dense cool gas 
(Section |g . 

2.2 Distance and size 

The small velocity vlsr of the shell and its longitude close to 
the Galactic centre make kinematic distance estimates uncer- 
tain as rotation curves here are poorly defined with significant 
velocity crowding. To deal with this uncertainty, we first cal- 
culate the kinematic distances and then constrain the results 
based on knowledge of supershell orig ins and dynami cs. We 
use a standard Galactic rotation curve (|Fich et alJll989h to es- 
timate the distance to the shell based on the systemic velocity 
of no 7 km s _1 , namely 

— = 1.00746 (^f) - 0.017112, (1) 
coo \ R J 

where uj is angular velocity, ujq is the angular velocity of the 
Sun, Ro is the Sun's distance from the Galactic centre and R is 

1 All velocities are quoted with respect to the Local Standard of 
Rest (LSR). 

2 We note that the systemic velocity here is taken to be the apparent 
central velocity of the shell, which also corresponds to the point of 
maximum absorption in its spectral profile (Section [3]) • 
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Figure 2. Hi velocity slices of GSH 006—15+7, over the velocity range 5—11 km s _1 in steps of ~ 1 km s -1 . The brightness ranges from 
10—80 K on a linear scale. Due to overlapping local emission at lower velocities, we see the structure of the shell most clearly at higher velocity 
although it is fainter. 
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Table 1. Basic parameters for GSH 006-15+7. 



Parameter 


Value 


Centre / 


6.0° 


Centre b 


-15.0° 


Central i; L sr 


+7 km s _1 


Distance 


1.5 ± 0.5 kpc 


Dimensions (angular) 


30° x 20° 


Dimensions (physical) 


790 x 520 pc 


^exp 


~ 8 km s _1 


E e 


~ 10 52 ergs 


Mass 


3 ± 2 x 10 6 M 


Age 


15 ± 5 Myr 



Galactocentric distance, using Ro = 8.3 kpc and Bp = 246 km 
s -1 in agreement with the most recent lite rature ( Ghez et al.l 
l2008l : lGillessen et alJl2009L iBovv et alJl2009h . For an input lon- 
gitude, latitude and velocity of / = 6°, b = —15°, v = 7 km s _1 
(the central coordinates of the shell and its systemic velocity), 
the resulting distance is d ~ 2.0 kpc for the near distance, and 
d rsj 15.5 kpc for the far distance. We do not further consider the 
far distance estimates due to the large angular size of the shell 
and corresponding physical size at such a distance. Due to the 
large angular scale of GSH 006-15+7, we concentrate on lon- 
gitude / = [7, 10]° and latitude b — [—2, 1]°. This is the region 
where we might expect to find a powering source based on shell 
morphology, and results in a near distance range of d — [1.2, 1.6] 
kpc. Rotation curve distance estimates generally provide a dis- 
tance that can be off by as much as 2— 3 kpc due to peculiar 
motions and local irregulariti es, an effect which is amplifie d 
particularly at low velocities ([Gomezl 2006; Ba ba et al1l2009h . 
If we consider the velocity gradient for the position of the shell 
over the velocity range [1, 10] km s _1 , this gradient equates 
to 0.3 kpc (km s -1 ) -1 , which indicates the magnitude of the 
error on the distance calculation. We hence use circumstantial 
factors to better constrain the distance to GSH 006—15+7. 

Looking towards the Galactic centre at longitudes between 
[0, 10]° on the near side of the Milky Way, we expect to cross 
the Sagittarius arm (1—2 kpc away), the Scutum-Centaurus 
arm (~ 3 kpc away), the Norma arm (> 4 kpc away) and 
the near 3 kpc arm (> 5 kpc away), where distances are esti- 
mated based on the best kn owledge of Milky Way s piral struc- 
ture using available data ([Churchwell et al . 2009). We note 
that there is a significant amount of star formation associated 
with a distance between 1.2 and 1.8 kpc, in the Sagittarius 
arm ([Steiman- Cameron et al.ll2010h , which agrees well with the 
kinematic distance of the shell. At a further distance of 3 kpc, 
if GSH 006—15+7 were located on the Scutum-Centaurus arm, 
we calculate the dimensions of the shell as ~ 1600 x 1000 pc. 
This would place the shell in the highest dimensions of known 
shells, especially for an object observed as a spatially coher- 
ent structure without clear signs of blowout. Based on this we 
believe that the shell is located near the Sagittarius arm and 
assign it a distance of 1.5 ± 0.5 kpc, with the errors reflecting 
the range of our kinematic distances. 

At an adopted distance of 1.5 kpc, the shell's angular di- 
mensions of 30° x 20° (latitude extent x longitude extent) 



(2) 



translate into physical dimensions of 790 x 520 pc. Based 
on these dimensions we can obtain a characteristic radius, 
~ 300+100 pc. We adopt this charac- 
teristic radius in order to estimate the properties of the super- 
shell and to compare it with the known supershell population. 
A characteristic radius of ~ 300 pc places the supershell in the 
largest 40 per cen t of supersh e lls, where the median shell ra - 
dius is - 210 pc (jHeiledl 19791 : iMcClure-Griffiths et alJ l2QQ2h . 
However, this sample is biased towards small nearby shells due 
to the difficulty in resolving small shells at any significant dis- 
tance. If we instead take both distance and radius into account 
by dividing the shell radius R s h by the shell distance <i s h, 

^sh 

where R s h and d s h are in parsecs, we see that the median 5 ~ 
0.05 and that GSH 006—15+7 lies within the top 3 per cent 
of supershells with 5 = 0.2, indicating that it is one of the 
nearest physically large supershells discovered so far. This is 
notable as the majority of other physically large supershells are 
quite distant in comparison with GSH 006—15+7, most likely 
due to the difficulty in identifying objects of large angular size. 
If the shell does extend nearly 800 pc off the Galactic plane 
as calculated, then we expect that it is almost reaching th e 
expected scale height of Galactic Hi ([Dickey &; Lockm an 1990). 



2.3 Mass, energy and age 

To estimate the mass of the shell, we first approximate its shape 
as a cone situated on top of a sphere. Summing the flux over 
only those velocity slices in which the shell is visible, we pro- 
duce a zeroth-moment map from which to extract measures of 
column density (assuming small optical depth), calculated as 



N H = 1.823 x 10 1 



7 



Tb dv cm 



(3) 



where Nn is the neutral atomic hydrogen column density, 
vo—vi represents its spectral extent and Tb is the observed 
brightness temperature. The mean column density was found 
by averaging the column density measured at different points 
along the wall of the shell. Assuming an axisymmetric wall 
width of the shell, we then use the background-corrected col- 
umn density of the shell. The background was approximated 
by sampling the column density away from the shell at several 
points and averaging to find the mean background, avoiding 
other regions of emission. This gives a mean corrected zeroth- 
moment of ~ 117 K km s _1 with a standard deviation of 
~ 15 K km s _1 . We convert the measured size of the cone 
and sphere from degrees to parsecs assuming a distance of 1.5 
kpc and derive a mean column density of Nh ~ 2 x 10 20 cm -2 . 

Once accounting for the presence of helium (a factor of 
1.4), our estimate for mass of the shell comes to M ~ 10 6 M© 
with a number density of ~ 1.3 atoms cm -3 in the shell and an 
ambient density of ~ 0.7 atoms cm -3 over the entire volume. 
Although there is foreground material present, using correla- 
tion between Hi and infrared data shows that the effect of this 
material on our mass estimate is negligible. We note that our 
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estimate here, typically used to estimate column densities of su- 
pershells, is affected by the necessary assumption that the gas 
is optically thin, and we correct for this assumption in Section 

El 

We derive an initial estimate of the total energy of the su- 
pershell by calculating the kinetic energy associated with its 
mass and velocity. We assume the mass to be as estimated 
(rsj 10 6 M©). Based on the large physical size of the shell and 
lack of a clear expansion signature in our data, we assume a 
late-evolution expansion velocity of 8 km s _1 , which is gener- 
ally assigned when shells are in their final phase of expansion 
and t he expansion velocity is comparable to random gas mo- 
tions (lHeiledll979h . This gives E K (0.5M^ xp ) - 6 x 10 50 ergs. 
If we assume that Er represents 20 per cen t of the total energy 
input due to efficiency o f energy transfer (Mc Crav &; Kafatod 
ll987l : lDawson et all 2008), then we would expect the total input 
energy to be ~ 3 x 10 51 ergs, which is likely to be an under- 
estimate of the total shell energy based on other supershells of 
comparable size. We can also obtain a theoretical estimate of 
the expansi on energy required to form the shell as describ ed an- 
alytically ([Weaver et al.lll977l : [McCrav fc Kafatod 1 19871 ). such 
that 



R = 97 pc(iY*/n 



1/5^3/5 



and 



5.7 km s-^iV./no) 175 ^ 275 



(4) 



(5) 



where R is the radius in parsecs, N* is the number of stars with 
M > 7 Mq (including both stars that have gone or will go su- 
pernova), no is the ambient density in cm -3 , and £7 is the age in 
units of 10 7 years. To incorporate the level of certainty to which 
we can specify no and R, we consider the range no = [0.5, 1.0] 
cm -3 (to account for uncertainties consistent with the ambient 
density predicted above) and R = [200, 400] pc (to account 
for the asymmetry of the shell). By solving Equations |4] and [5] 
simultaneously, we find that £7 = [15, 30] Myr and N* — [6, 
47]. This corresponds to an input energy of [1, 10] xlO 52 ergs, 
assuming each star contributes 10 51 ergs from stellar winds and 
10 51 ergs from going supernova. Overall these energy estimates 
suggest a formation energy of the order of 10 52 ergs, which is 
most likely the degree of accuracy to which we can predict the 
formation energy without further data or modelling. 

For comparison to the overall distributi on of super- 
shell energy as performed above for radius ([Heiled Il979l ; 
iMcClure- Griffiths et al. 2002), we find a median known super- 
shell energy of ~ 9x 10 51 ergs which places GSH 006—15+7 very 
close to the average formation energy of known supershells, in 
rough agreement with the percentile of the characteristic ra- 
dius. Estimating the age of a shell is very difficult in the ab- 
sence of a current powering source, and so the age represented 
by £7 is highly uncertain. We can however assume an upper age 
limit of around 20 Myr, after which the combination of ambient 
pressure at high scale height, galactic shear and gravitational 
deceleration likely destroys mo st coherent structure in a su- 
pershell and leads to breakou t (|Tenorio-Tagle &; Bodenheimerl 
1988: lMashchenko et al.lll999h . 




GLON (degrees) 

Figure 3. Map of Ha emission (greyscale) near GSH 006—15+7 
(Finkbeiner 2003). In all images, the contours show Hi emission (yel- 
low, at 25, 30, 35 K) at ~ 7 km s _1 (top) and ~ 10 km s _1 (bottom), 
100 /im dust (red, at 8, 9, 11 MJy/Sr) and Ha emission (black, at 3, 
6, 9 R). At 7 km s _1 , we see strong Hl/IR correlation. 



2.4 Multiwavelength properties 

We investigated available multiwavelength data to determine if 
other emission processes are evident. Previously detected emis- 
sion in supershells include Ha, s oft X-rays, polarised radio con- 
tinuum an d 100 /im emission (Heil es et al.l Il999l : IWest et al. 
l2007l : lMcClure-Griffiths et al.ll200lh. We find no convincing ev- 
idence of soft X-rays (|Snowden et al.lll995h (due to co nfusion) 
or 1420 MHz polarised emission (|Testori et al.l l2Q08h associ- 
ated with the shell. We also looked at the Ha data from the 
South ern H- Alpha Sky Survey Atlas ( SHASSA;lGaustad et al.l 
l200lh Ha imag ing survey co mpiled bv lFinkbeinerl (|2003T ) . Based 
on extinction flS chlegel et al. 1998) and optical depth estimates 
(|Madsen &; Revnolddl2005r )in this region of the Galaxy, we can 
place an upper limit on the intensity of any possible Ha emis- 
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Figure 4. Ha spectrum obtained with WHAM. This is an on— off 
spectrum, as described in the text. An Hi spectrum from an aver- 
age of the GASS beams falling within the WHAM beam is shown 
with the dashed red line; an average of Hi spectra from GASS beams 
falling within a beam the size of the WHAM beam but centred in a 
nearby knot of bright Hi is shown with the dotted red line. The geo- 
coronal Ha line is incompletely subtracted near i;lsr = —20 km s -1 . 



sion by noting that we see no emission within the main shell 
structure (b ~ —20°) or along the shell walls (b ~ —25°) at lim- 
its of [1.2±0.5, 2.0±0.5] R0 respectively, which, given estimated 
optical depths of r = [0.23, 0.46] at 656.26 nm, corresponds to 
maximum surface brightnesses of [1.5 =b 0.6, 3.1 =b 0.8] R. 

At higher latitudes where the shell may be fragmenting, 
we note an evident Ha structure in the SHASSA data at 
/ ~ 8°, 6 ~ —30° with an angular extent of ~ 2° and a peak 
brightness of ^ 15 R. This structure is on the interior of the 
fragmentary Hi emission seen at high latitudes (described in 
Section I2.1|) . We find that the Hi at lower velocities correlates 
well w ith IRIS 100 /xm emission ([Miville-Deschenes &; Lagachd 
2005). Fig. [3] shows the emission in Hi, Ha and 100 fim at 
different velocities, showing correlation between the shell and 
the Ha feature. To verify that the morphological association 
of the Ha and Hi corresponds to physical association, we ob- 
tained an Ha spectrum with the Wisc onsin H- Alpha Map- 
per (WHAM; lHaffner et ail l200l l2010h towards the region 
of the cloud with the brightest Ha emission in the SHASSA 
data, (l,b) — (8.70,-29.85). The data were obtained from 
Cerro Tololo Interamerican Observatory on the 16th Novem- 
ber 2011. WHAM integrates all emission within its l°-diameter 
beam. We also obtained an off-source Ha spectrum centred at 
(7.03, —33.92). We fit a single Gaussian modelling the geocoro- 
nal Ha line to both the on- and off-source spectra and then 
subtracted the geocorona-subtracted off spectrum from the 
geocorona-subtracted on spectrum to remove the faint, mostly- 
unidentified atmospheric lines that dominat e the baseline at th e 
surface brightness sensitivity of WHAM (|Hausen et al.ll2002i ). 
We also used the geocoronal line to calibra te the velocity scal e 
to an accuracy of less than 0.5 km s _1 ([Haffner et al.l 120031 ). 
The resulting spectrum consists of entirely astronomical emis- 



3 One Rayleigh (R) = 10 6 /47r photons s 1 cm 



sion and is shown in Fig. 2] The emission is well-fit by a single 
Gaussian component with a mean LSR velocity of 6.7 km s _1 , 
a full width at half-maximum of 27 km s" 1 , and an inten- 
sity of 5 R. Although the relatively large width of the Ha 
line precludes definitive separation from local gas, the single- 
component Gaussian profile and the good agreement of the 
velocity centroid of the Ha emission with the Hi data strongly 
support the association between the Ha and the Hi-emitting 
gas which is suggested by the morphology from the SHASSA 
image. 

While Ha is observed in nearly all directions at high lat- 
itude (jHaffner et al.l 12003). multiple emission mechanisms are 
likely responsible in different directions. Ha well off the Galac- 
tic disk is thought to be emitted either by gas which is pho- 
toionised in situ by ionising phot ons from OB stars in the disk 
which have escaped th e plane ([Reynolds! Il990l ; lHaffner et al.l 
2009: IWood et alJfeOld ) or to be Ha emission from closer to the 
lane which has been scattered off o f high-latitude dust clouds 
IWitt et al.1 l20ld : ISeon et al.l l201lh . Ha can also be emitted 
shock-ionised gas or by gas photoionised by a harder (soft X- 
ray) radiation field such as might be produced by million-degree 
gas in the interior of the supershell. With future observations 
of the [O in] A5007 and [S n] A6716 e mission lines in combi- 
nation with photoio nisation modelling ([Wood h, Mathis 2004; 
iMadsen et al.l l2006), we may be able to discriminate between 
emission mechanisms in this Ha feature and determine how the 
Ha relates to the supershell. 



3 HI SELF-ABSORPTION IN GSH 006-15+7 

A unique feature of GSH 006—15+7 is that we see the up- 
per right shell wall transition clearly from Hi emission into 
absorption at b > —6.5°, with the deepest absorption occur- 
ring at the systemic velocity v ~ 7 km s _1 . Fig. [5] shows the 
transitioning wall at this systemic velocity, traced with a line. 
The region where we might expect to find the transition point 
is marked with a cross, showing a shift from emission to ab- 
sorpti on. This phenomenon is known as H i self-absorption, or 
HiSA (lHeeschenlll955l : lGibson et alJl200Ch . An absorption pro- 
file (such as the example spectrum shown in Fig. [6} can in- 
dicate the presence of cool atomic gas. We expect to see Hi 
absorption against a diffuse background Tb g that has a higher 
brightness temperature than the spin temperature T s of the 
foreground gas. However, the detection of absorption is not 
generally enough to constrain the spin (or excitation) tempera- 
ture of the gas and so it is usually assumed, constrained based 
on the presence of other m olecules or estimated via line width 
fitting ([Gibson et al]l2000h . In special cases, a transition from 
Hi emission to absor ption allows T s to be measured directly at 
the transition point (|Kertonll2005h . 

We apply HiSA diagnostics to this feature below in order 
to estimate spin temperature T s , optical depth r and column 
density Nh in the shell. We note that there is some spatial 
proximity of GSH 006—15+7 to the local self- absorbing Riegel- 
Crutch er cloud at a distance of ' y 100 pc towards the Galactic 
centre ([McClure-Griffiths et alj |2006). with their angular sep- 
aration varying from a few to about 15 degrees. Taking the 
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Figure 5. Image of the region of the shell used for HlSA analysis, 
at the systemic velocity of ~ 7 km s -1 . We scale the data linearly 
from 20 K to 70 K. The line traces the position of the shell wall and 
the cross marks the region where GSH 006—15+7 transitions from 
absorption to emission. 



point of maximum absorption to represent the systemic ve- 
locity of Riegel-Crutcher gives a peak absorption velocity of 
^5.8 km s _1 , compared to the peak absorption velocity of 
GSH 006-15+7 at - 7.2 km s" 1 . The two objects also show 
differing absorption profiles and are morphologically distinct, 
but with some evident overlapping absorption at latitudes > 0°. 
In our HlSA analysis, we use only latitudes < 0° in order to 
avoid any confusion. 

In a simplified case of absorption in the context of radiative 
transfer, we can define an 'on-source' brightness temperature 
(Tb) where we observe the absorption through the absorbing 
cloud (also known as the observed temperature, which in this 
case is the continuum-subtracted Hi brightness temperature), 
and an 'off-source' brightness temperature (Tb g ) where we es- 
timate the background without the absorption present. Under 
the limiting, but reasonable, assumptions that an y background 
or fo reground gas has negligible optical depth ([Gibson et al.l 
2000), the difference between the on and off is 



-T hs = (T B -T c -pT bs )(l-e- T ) 



(6) 



where T s is spin temperature of the absorbing cloud, p is the 
fraction of Hi emission originating behind the HlSA cloud (as 
opposed to i n front of th e cloud), and T c is the mean continuum 
background (|Feldtlll993l ). A value of p — 1 assumes that all Hi 
emission originates from behind an object seen in absorption. 
Although we know there is some foreground material towards 
GSH 006—15+7, setting p = 1 is a reasonable assumption given 
the relatively small path length to the shell. In the case of 
a transition from emission to absorption, the transition point 
marks the place where the background Tb g and the source Tb 
are equal in magnitude; in this case Equation [6] simplifies to 



Figure 6. Example spectrum at a single position (I = 5.395, b = 
—3.236) showing the absorption due to GSH 006—15+7. The maxi- 
mum absorption is predominantly seen at the systemic velocity of 7 
km s _1 along the wall of the shell. The maxima/minima (circles) and 
estimated FWHM (squares) are shown overplotted on the spectrum 
and are used in Section 13.11 to estimate a line- width of ~ 2 km s _1 
corresponding to a spin temperature of ~ 40 K. 



Equation thus allows direct calculation of the spin temper- 
ature under certain circumstances at the transition point in 
GSH 006—15+7, giving insight into the local properties of the 
Hi gas. Once we estimate the spin temperature T s , it is possi- 
ble to substitute the values for T&, Tb g , T s and T c to solve for 
optical depth, using 



r = -In 1 



tbg 



(8) 



Using Equation [8] along the wall of GSH 006—15+7 under the 
assumption of constant T s , we can calculate r and use this to 
estimate the corresponding column density, using 



1.823 x 10^7; / rdv 



(9) 



which, assuming 
plifies to 



t T s Av 
C ' 



, Gaussian form for the integral of r(v), sim- 



(10) 



T e =pT bs + T c . 



(7) 



where Nn is the column density, To is the line centre opacity, 
T s is the spin temperature, A v is the line profile width and C 
is a standard correction term dGibson et al.ll20Q0h [5.2 x 10" 19 
cm 2 K km s _1 ] (|Dickev &; Lockmanf 1990l ) From here, we can 
redo our previous mass calculation under the assumption that 
this newly estimated column density better represents the true 
value as the unlikely assumption of thin optical depth is no 
longer made. 
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Spectral HISA (step = 0.2°) 



on source 
off source 
continuum 



Position Along Feature (deg) 

Figure 7. Spectrally-derived variation in brightness temperature 
on/off the shell arm. This plots the maximum absorption at each 
latitude (diamonds) versus the background temperature (squares) 
and the continuum background (circles, Calabretta et al, in prepa- 
ration), with the amount of absorption maximised along the wall of 
the shell. The horizontal axis is zeroed at b = —2°. The grey-shaded 
region represents the extent to which we believe the background is 
uncertain (~ 10 K), which is incorporated into the calculations. The 
transition point occurs where the curves meet, but the consequent 
emission of the shell is not clear due to the difficulty of spectrally 
separating emission due to GSH 006—15+7 and other emission. 



3.1 HiSA results for GSH 006-15+7 

We examine the absorption to emission transition in velocity- 
longitude space along the wall of the shell at 0.2° increments 
over the range b = [—7.0, —2.0]°, obtaining the spectrum at the 
velocity of maximum absorption for each latitude point. We re- 
sample each spectrum in order to ease the location of minima 
and maxima along the curve, taking into account the velocity 
resolution of ^ 1 km s _1 in estimating the velocity width of 
each profile. By obtaining the minimum and maximum points 
that define the absorption, we adopt these as our on-source 
(Tb) and off-source (Tb g ) values. We determine the continuum 
background (T c ) from Hi Parkes All Sky Survey (HiPASS) con- 
tinuum data at 1.4 GHz corrected for extended emission (Cal- 
abretta et al, in preparation). The results are plotted in Fig. [3 
showing the on-source Tb, the averaged off-source Tb g and the 
background continuum T c . 

We assume the simplest case for the background by taking 
the average of each pair of maxima to represent the amount 
to which the emission is absorbed. To investigate the error on 
our assumed background, we examined longitude profiles at 
7 km s _1 at the latitude of each point used in our analysis 
and find the peak brightness temperature on the right side of 
the spatially-apparent absorption (outside the shell wall). This 
allows us to check how well the background derived from the 
velocity spectra agrees with a spatially-determined background 
temperature. We find a positive trend indicating that the spec- 
tral profile does underestimate the background, with a peak 



deviation of ~ 10 K and a mean deviation of ~ 4 K. We thus 
assume that our background estimate could be as much as ~ 10 
K underestimated and factor this into our predictions below. 

The spectral method we have used, although reliable, does 
not easily allow visualisation of the emission of GSH 006—15+7 
after the transition point because it is difficult to separate the 
shell emission from background emission without stable back- 
ground spectral profiles that can be subtracted. Due to this, 
we cannot use Equation [Tito tightly constrain the spin temper- 
ature at the transition point, but we can set an upper limit to 
the spin temperature by considering the combined Tb and T c 
values along the wall, 



T b + T c; 



(11) 



as in our simple radiative transfer this will always be an over- 
estimate as there are likely to be sources contributing to the 
emission along the line of sight that are not taken into account. 
We thus calculate T s?ma x systematically along the shell wall and 
then average the results to give a mean maximum spin temper- 
ature of Ts^max ~ 40 K. Any fluctuations in spin temperature 
due to the local properties of the gas in the shell wall are likely 
to be moderated by using the mean value. To check the value 
of T s , we estimate the width of the absorption profile in Fig. 
[6] and find the FWHM - 2 km s" 1 . In this case, the resulting 
thermal velocity dispersion a v is FWHM/2.36 ~ 0.8 km s _1 
with a kinetic temperature of 



mu<Jy 
k ' 



(12) 



where mu is the mass of hydrogen and k is Boltzmann's con- 
stant. This yields T k - 90 K. For temperatures T k < 1000 K in 
the cold neutral medium (CN M), T k is essentially equal to T s 
(|Lisztll200ll ; iRov et al.l 12006). However, spectral modelling of 
Hi absorption and emission profiles demonstrates that signifi- 
cant line broadening relative to the spin temperature at peak 
opacity r m ax is generally seen (|Heiles &; Trolandll2003l ) and can 
indicate either or both of significant internal temperature struc- 
ture in the gas and a turbulent contribution to the line-width 
measurement. Our spin temperature estimate of 40 K implies 
excess line-broadening of a v ~ 0.6 km s _1 , which is fully con- 
sistent with values typically found in the ISM. 

Using Equation [8] along the absorbing wall of 
GSH 006—15+7 and a spin temperature of 40 K gives 
values of r in the range [1.8, 5.0] with a mean of r ~ 3.1. The 
potential underestimation of the background of up to ~ 10 K 
is used to calculate the upper limit in the given ranges. We 
use the optical depths to estimate the corresponding column 
density, which ranges over [2.7, 7.7] xlO 20 cm -2 with a mean 
value of Nu ~ 4.7 x 10 20 cm -2 . We compare this to the 
value estimated from the zeroth-moment map in Section 12.31 
of Nu ^ 2 x 10 20 cm -2 , and note that the mean column 
density is more than double the previous estimate, which is to 
be expected due to eliminating the assumption of negligible 
optical depth. 

If we assume that our new estimate better reflects the 
true column density and substitute this into the calculation for 
mass, we find that M ranges over [1.8, 5.1] x 10 6 Mq with a cor- 
responding number density of [1.7, 4.8] atoms cm -3 . This is in 
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comparison to our previously determined estimate of ~ 10 6 Mq 
and ^1.3 atom cm -3 , where the assumption of zero optical 
depth leads to a significant underestimation of the shell mass. 
We note that the application of a single opacity correction fac- 
tor to the entire supershell may be erroneous as it would be 
plausible to have systematically lower densities and opacities 
at larger off-plane distances, while the opacity measurement is 
done closer to the plane where these would be maximised. To 
counter the effect of this, we use the range of resultant masses 
from using optically thin versus HISA as a representation of 
the error on the mass such that M = 3 =b 2 x 10 6 M©. 



4 EVOLUTION OF GSH 006-15+7 

To investigate further the nature of the shell, we have searched 
for potential stellar powering sources and applied simple mod- 
elling in order to constrain the origin of GSH 006—15+7. This 
involves estimating the amount of energy deposited by candi- 
date sources through stellar winds and supernovae in order to 
assess whether these regions could form GSH 006—15+7 based 
on the energy predictions we made in Section 12.31 where we 
found N* ~ 26 and E ~ 10 52 ergs. 



4.1 Energy estimation and supernovae 

We use mult i wavelength data and published databases to in- 
vestigate the region of the Galactic plane where the walls of 
GSH 006—15+7 meet to search for poweri ng sources. Of the 
nearby OB associations ([Humphreys! 1 19781 ), Sagittarius (Sgr) 
OB 1 is a likely candidate based on coordinates as well as 
its reported distance of 1.58 kpc. The 13 high-mass stars of 
Sgr OB 1 are spread over 10° of longitude, which makes it dif- 
ficult to assess their true association. At the reported distance 
of rsj 1.6 kpc, the mean longitudinal drift of each star from 
the centre of the Sgr OB 1 association corresponds to roughly 
40 pc, in agre ement with an average stellar drift velocity of 
- 5 km s" 1 (|McCrav fe Kafatod [l987h and an age of around 
10-20 Myr for GSH 006-15+7. We show their distribution 
in Fig. [8l as well as the distribution of the nearby and possi- 
bly associated open clusters NGC 6514, NGC 6530, and NGC 
6531 (with which some of the stars of Sgr OB 1 are associ- 
at ed) . The general area has also b een referred to as Sagittarius 
I (Sta lbovskv fe Shevchenkolll98lh . incorporating a number of 
clusters and nebulae. Overall, the entire region shows evidence 
of high star formation activity, and so we base our initial en- 
ergy calculations on the stellar winds of the known high-mass 
stars. 

We calculate stellar wind estimates using Sgr OB 1 , 
of which stellar types ar e well known ([Mason et al ] Il998l : 
iMafz-Apellaniz et alj|2004 ) and can be used to determine ef- 
fective temperature, luminosity, mass, maximum age, mass loss 
rate, wind terminal velocity and wind luminosity. We estimate 
the stellar wind mechanical luminosity as Lsw = |MV^, based 
on the observed spectral type which allows derivation of mass 
loss rate M and wind terminal velocity Voo • We match t he stel- 
lar ty pes earlier than B0 to those for the Carina Nebula ( Smith! 



generally not included in models) by raising them to the next 
highest class (i.e. I and III), which provides an overestimate for 
their energy contribution. We have also assumed that the age 
of each star is the lifetime of a star with its estimated mass, 
which gives us an upper limit. Given the resulting wind lumi- 
nosity, it is then possible to estimate total energy input from 
these stars due to their stellar winds. 

It is expected in studies of large old shells that many of 
their forming high-mass stars will have exploded as supernovae 
(SNe), as their lifetimes are on the order of a few to several 
million years ([Chiosi et al.l ll992). As such the energy calcu- 
lated from stellar winds alone of any powering region is likely 
to underestimate its true energy contribution, and so we assess 
the contribution of stars that have already gone supernovaEl 
We can estimate the number of supernovae based on the ini- 
tial mass function (IMF), a measure of mass distribution. This 
can be readily done in the case of stellar clusters with a known 
IMF, as opposed to OB associations whose stars can be scat- 
tered across large regions of space and possibly even associated 
with different stellar clusters (as in the case of Sgr OB 1). We 
perform a cursory IMF calculation for Sgr OB 1 but emphasise 
that there are strong assumptions impacting this estimate. For 
this reason, we focus the IMF calculation on the open clusters 
NGC 6514, NGC 6530 and NGC 6531 due to their proximity to 
GSH 006—15+7, distance close to 1.5 kpc and stellar content, 
noting that several of the stars of Sgr OB 1 are associated with 
these clusters. We do not include the other posit ionally nearby 
clusters NGC 6546 (due to its large estimated age exceeding 
the shell's age) or Bochum 14 (due to its lack of high-mass 
stars and poorly known stellar composition). 

For the clusters, we use UBV photometry and stellar clas- 
sification of members to calculate effective temperature (log 

of each star (Mass ev et alj 119951 ), 
calculate the absolute magnitude (My) 



T e ff), bolometri c correction (BC) and V-band extinction (Ay ) 

using these values to then 



bolo metric magnitud e 
(Mboi) and hence luminosity (log L/L©) ([Crowtherl I 12004 ). 
Masses are estimated using this resulting lu minosity via the 
mass- luminosity relationship ( E ddingtonlll924 ) , with awareness 
that this relation is generally valid only for main sequence stars. 
We make the assumption that the photometry is complete from 
the minimum high-mass star calculated (M m i n ) to the maxi- 
mum high-mass star calculated (M max ), with any stars higher 
than this mass already gone sup ernova. We a dopt a minimum 
possible stellar mass of 0.5 M© (Kroupa 200l|) and a maximum 
possible stellar mass of 150 M (Figer 2005), and integrate over 
the range assumed complete 



F N (M min < M < M max ) 



M 1T 



M Q dM 



Jo 5 M * dM 



(13) 



where Fn is the number fraction of stars in the given range 
with respect to the total and a is the slope of the initial mass 
function. We can therefore estimate the total number of cluster 
stars by dividing the number of stars N in this range by Fn- 
Then, we estimate the number fraction of exploded stars in the 
cluster, 



lar ty 
2006) 



and account for luminosity classes II and IV (which are 4 We refer to these stars as exploded stars. 
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F N (M > M n 



M a dM 



(14) 



and multiply by the total number of stars to obtain an estimate 
of the total number of exploded stars wh i ch we attribute to 
super novae (jBruhweiler et al.lll98Ql : lHeiled ll987). We assume 
that each exploded star has contributed 10 51 ergs of energy 
through stellar winds over its lifetime and 10 51 ergs via its 
supernova, giving an energy contribution of 2 x 10 51 ergs per 
exploded star. 



4.2 Cluster parameters 

We obtain UBV photometry and stellar types for each clus- 
ter from the online WEBDA (Web Open Cluster Database 
(jMermilliod fe Paunzenl [2003) . As we are necessarily assum- 
ing that these clusters contributed to the formation of 
GSH 006—15+7, we adopt a distance of 1.5 kpc to all clusters. 
We ta ke the age estimates from the catalogue of Dia s et all 
(2002), and obtain values of mean reddening from the litera- 
ture. We obtain values for the IMF slope (a) of each cluster 
where available, ot herwise we ado pt the traditional Salpeter 
value of a — — 2.35 (|Salpeterlll955h . We use in our calculations 
a total of 60 stars from the three clusters described below. 



4.2.1 NGC 6514 

This object is also known as the Trifid Nebula or M20, housing 
an open cluster of around 360 known stars as well as various 
types of nebulae and dense regions of dust. The distances re- 
ported for this object are quite varied, from 0.8 kpc to over 2 
kpc, however we adopt a distance of 1.5 kpc, the Dias cata- 
logue a^e_o£^3 : 3_^^rarid a mean reddening of E(B — V) = 
0.23 (jOgura fe Ishidall 19751 ). We use an IMF slope of a = -2.35 
for this cluster. 



4.2.2 NGC 6530 

NGC 6530 is a well-studied young, active star-forming cluster 
associated with the Lagoon Nebu la, described as the core clus- 
ter of the Sgr OB 1 association (jChen et alJl2Q07h and housing 
roughly 350 known stars. Of the three clusters, it is the most 
active and also houses the highest number of O stars. Based 
on th e literature, we a dopt a mean reddening of E(B — V) = 
0.35 (Sung et al. 2000) for this cluster, the Dias catalog ue age 
of 7.4 Myr and an IMF slope of a = -2.22 (jPrisinzano et al.1 
l2005h . 



4.2.3 NGC 6531 

This cluster is not particularly well studied due to its small 
number of detected stars (around 100), lack of numerous high- 
mass stars, nebulae or bright st ar forming regions. It has re - 
ported distance mo duli of 11.35 (Mermilli od &; Maederl [l986h , 
10.70 ([Forbeslll996l). 10.5 (IPark et al.ll2Q0lh and most recently 
10.47 (Mc Swain &; Giesl l2005), giving a distance range of 1.2 to 
1.9 kpc with a mean of 1.4 kpc. Its age has been estimated with 
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Figure 8. The assumed source region of GSH 006—15+7, showing 
the stars of Sagittarius OB 1 (red crosses) and the open clusters (red 
boxes) used for stellar wind estimates and for estimating the energy 
contribution of supernovae overlaid on a Hi velocity slice at ~ 5 km 
s _1 (where the cavity these objects occupy is most clearly visible). 
The black contour is at 55 K to show the outline of the shell and the 
scale is linear from 30 K to 70 K. 



a reported range of 8—11 Myr. We adopt the Dias catalogue 
age of 11.7 Myr, a mean reddening of E (B — V) — 0.28 for this 
object and an IMF slope of a = —2.29 



4.3 Modelling the energy of the shell: stellar winds 

In the case of Sgr OB 1, we estimate ~3x 10 51 ergs of stellar 
wind energy, too low to account solely for the formation of 
GSH 006—15+7. Similarly we obtain a total cluster stellar wind 
energy of ~ 1 x 10 51 ergs, which we expect because the clusters 
have fewer high- mass stars overall. As only 20 per cent of these 
winds would have contributed to the expansion of the shell 
and our energy prediction was 10 52 ergs, we conclude that the 
stellar winds of Sgr OB 1 or the clusters could not be the sole 
powering source of GSH 006—15+7, which is not necessarily 
surprising given the size and age of the shell. We therefore 
assume that the majority of the expansion energy must have 
come from supernovae that are now invisible. 



4.4 Modelling the energy of the shell: supernovae 

We now estimate the amount of energy we could derive based 
on exploded stars from the clusters we have examined. As be- 
fore we adopt a dista nce of 1.5 kpc, v alues of reddening from 
the literature and the lDias et al.l (|2002h catalogue cluster age 
estimate. The results are summarised in Table [2] We obtain 
a maximum energy input from supernovae and stellar winds 
(from exploded stars) of around ~ 2 x 10 52 ergs total. But us- 
ing the same conversion efficiency of supernova energy used in 
Section [231 of 20 per cent gives us only EsNe ~5x 10 51 ergs. 

For comparison, we perform the same calculation for the 
stars of Sgr OB 1 (for which we have more reliable masses) 
but emphasise that this approach implicitly assumes that the 
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Cluster C 1 ) 


M > 8 M ( 2 ) 


M min (M ) ( 3 ) 


M max (M ) ( 4 ) 


a (5) 


Cluster stars ( 6 ) 


# SNe ( ? ) 


Energy ( 8 ) 




NGC 6514 


[4, 11] 


[5, 7] 


[23, 65] 


-2.35 


[360, 520] 


[2, 1] 


[4.0, 2.0] xlO 51 


ergs 


NGC 6530 


[20, 30] 


[5, 8] 


[28, 83] 


-2.22 


[590, 970] 


[4, 1] 


[8.0, 2.0] xlO 51 


ergs 


NGC 6531 


[5, 9] 


[4, 6] 


[17, 42] 


-2.29 


[260, 400] 


[3, 2] 


[6.0, 4.0] xlO 51 


ergs 


Total 












[9, 4] 


[1.8, 0.8] xlO 52 


ergs 


Sgr OB 1 


[13, 13] 


[17, 19] 


[43, 56] 


-2.35 


[2100, 2300] 


[4, 3] 


[1.0, 0.6] xlO 52 


ergs 



Table 2. Parameters of each cluster used to model IMF and resulting values. The columns are: (1) cluster name, (2) number of stars greater 
than 8 solar masses, (3) minimum mass used, (4) maximum mass used, (5) initial mass function slope, (6) extrapolated number of cluster stars 
based on IMF predicted number fraction, (7) number of supernovae estimated based on the mass distribution and (8) corresponding energy of 
these stars, assuming 10 51 ergs from stellar winds and 10 51 ergs from the supernova. In all cases for the clusters, the first number represents a 
mass po wer of 4, and the second number repr esents a mass power of 3. In the case of Sgr OB 1, the numbers represent masses from different 
sources ([Martins et al.l l2005: Wei dner fe Vinkll2010h . 



association stars formed from the same mass distribution. This 
is likely to be an o ver-simplification. We adopt mass estimates 
from two sources ([Martins et al J 120051 : IWeidner fc Vinkll2Q10h 
and use this to give a likely energy range. We use the standard 
Salpeter value of a = —2.35 for the IMF slope and assume that 
the association stars are complete from the minimum detected 
mass (17 M©/19 M©) to the maximum detected mass (43 
M /56 M©), where the two numbers represent mass estimates 
from different sources (Martins et al . 2005; Weidn er fc Vinkl 
2010). This results in an energy range of 0.6—1.0 xlO ergs 
total, which is essentially in agreement with the cluster es- 
timate as suggesting a maximum converted energy on order 
- 2 x 10 51 ergs. 

We find here that the formation of GSH 006—15+7 cannot 
be accounted for solely by the supernovae and stellar winds of 
exploded stars from either the clusters or Sgr OB 1. The com- 
bination of stellar wind energy found in Section 14.31 (^ w ind ~ 
5 x 10 50 ergs) with that of the exploded stars (EsNe ~ 5 x 10 51 
ergs) is close to 10 52 ergs (with the dominant contribution being 
from exploded stars). We also expect that our rough assump- 
tion of 2 x 10 51 ergs per exploded star may in fact be an under- 
estimate of the true energy contribution of the very high-mass 
stars, since the winds are likely to e xceed the supernova energy 
for high-mass stars earlier than 06 (Abbott 1982). It therefore 
seems reasonable to assume that the shell was formed by the 
stars in and or around the Sgr OB 1 association. 



5 CONCLUSION 

We have reported on the discovery, properties and analysis of a 
new Milky Way supershell GSH 006-15+7. This object spans 
longitudes of / = [356, 16]° and latitudes of b = [—28, 2]°, which 
translates into physical dimensions of (790 ± 260) x (520 =b 
175) pc at an estimated distance of 1.5 ± 0.5 kpc. The shell 
is elongated in the latitude direction, which is possibly due to 
a density gradient away from the Galactic plane. We estimate 
the dynamical age to be between 7—15 Myr, with the coherent 
structure of the shell suggesting an age less than 20 Myr based 
on known supershells in the Galaxy. Energy estimates suggest a 
formation energy most likely ~ 10 52 ergs, with a mass estimate, 
corrected for optical depth, of 3 ± 2 x 10 6 M©. 



We have found no convincing evidence of an upper shell 
counterpart. In the case of an absent upper shell complementing 
GSH 006—15+7, we can interpret this as the result of an asym- 
metric outflow, due possibly to inhomogeneous gas surround- 
ing the source or a source region slightly below the Galactic 
plane which preferentially outflows in the direction of the low- 
est density. An offset of even 100 pc from the Galactic plane, 
depending on the vertical density distribution, is known to have 
a significant effect o n the direction of outflow an d shape of the 
resulting supershell ([Chevalier &; Gardner! 1 1974[ ). 

We have investigated a self- absorbing feature along the 
shell wall in order to place a more realistic constraint on our 
estimates of shell mass as well as giving an indication of gas 
temperature. We have found that the gas in the supershell is 
comparatively cool with a spin temperature of ~ 40 K as well 
as being quite optically thick (and hence dense) with r ~ 3 
which is supportive of recent studies showing correlation be- 
tween supershel ls and molecular ga s formation, which requires 
dense cool gas (|Dawson et al.l 120081 ; iFukui &; Kawamural l2010l ; 
iDawson et alJfeOllh . 

We see strong evidence of past star formation in the re- 
gion where we would expect to find a powering source based 
on the morphology and physical properties of the supershell. 
We also find some evidence in favour of a powering source 
for GSH 006—15+7 via the OB association Sgr OB 1 and the 
nearby open clusters NGC 6514, 6530 and 6531. In analysing 
the energy contribution from stellar winds from existing stars 
and the winds/supernovae of exploded stars from Sgr OB 1 
and the clusters, the estimated combined ~ 5 x 10 51 ergs en- 
ergy is close to the estimated requirement of ~ 10 52 ergs of 
formation energy, and is within the errors in our analysis. The 
absence of presently observable high-mass stars contributing 
significant stellar winds suggests an expansion partially, if not 
mostly, driven by the winds and supernovae of exploded stars, 
which likely originated as part of the same region occupied by 
Sgr OB 1 and the clusters. 

Based on our analysis, we suggest that GSH 006—15+7 is 
in the transition stage between supershell and chimney. The 
filamentary structures we see in the morphology of the shell at 
high latitudes, evident in Hi and IR with possible associated Ha 
emission, may be an indication that fragmentation and blowout 
is already beginning to take place. 
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